Local correlations of different eigenfunctions in a disordered wire 
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We calculate the correlator of the local density of states {pe{ri)ps+u:{r2)) in quasi-one-dimensional 
disordered wires in a magnetic field, assuming that |ri — r2| is much smaller than the localization 
length. This amounts to finding the zero mode of the transfer-matrix Hamiltonian for the su- 
persymmetric cr-model, which is done exactly by the mapping to the three-dimensional Coulomb 
problem. Both the regimes of level repulsion and level attraction are obtained, depending on |ri— r2|. 
We demonstrate that the correlations of different eigenfunctions in the quasi-one-dimensional and 
strictly one-dimensional cases are dissimilar. 
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cd 



I 

O 

o 



> 

m 

o 
o 



I 

O 

o 



% 



Modern physics knows a number of paradigmatic mod- 
els which describe properties of various, seemingly unre- 
lated, phenomena. Anderson localization in disordered 
systems with one-dimensional geometry is one of such 
models. Historically, it played an important role in set- 
ting down the concept of quantum localization. The 
model of a one-dimensional (ID) chain was the first ex- 
ample of the system where localization takes place at 
arbitrary weak disorder [l|, |2|. Later on Thouless p] 
argued that all the states are localized also in disor- 
dered quasi- one- dimensional (QID) systems, supporting 
a large number iV ^ 1 of propagating transverse modes. 

During the last three decades it had been realized that 
many problems of condensed matter physics and quan- 
tum chaos can be mapped onto QID localization. The 
most natural example is particle propagation in a disor- 
dered wire, both in the limit of weak jj| and strong (gran- 
ular) 5] disorder. Besides that, the problem of random 
banded matrices with a large bandwidth can be mapped 
[6j onto the same model. The problem of the quantum 
(5-kicked rotor '7\ 8] whose evolution operator looks like 
a quasi-random banded matrix also belongs [9|, llCl t o the 
QID universality class. Recently it was argued [ll| that 
dynamic localization in quantum dots which can be de- 
scribed by time-dependent random matrices is equivalent 
to QID localization as well. 

Despite the fact that mathematical description of 
strictly ID disordered systems is rather involved ll2Lll,ll |. 
our knowledge about them is quite complete. In this case, 
the localization length ^id coincides with the mean free 
path I. The low- frequency conductivity follows the Mott- 
Berezinsky law: Kea{uj) ex cj^ In (l/wr) |l2l. Il4l|. where 
T is the elastic time. The correlation function of the lo- 
cal density of states (LDOS) was considered by Gor'kov, 
Dorokhov and Prigara ,15] who have shown that in the 
limit of small energy separation, ujt ^ 1, the eigenstates 
are uncorrelated at r ^ ^id ln(l/a;r), exhibit nearly per- 
fect level repulsion at $id <C r ^ ,^id ln(l/wT), with a 
tendency to level attraction at r <C i^id (for a discussion 
of these results using Mott's arguments, see Ref. I1(t|). 



In the QID case, the localization length ^ ^ Nl 17] is 
much larger than the mean free path, that allows to for- 
mulate the problem of QID localization on the language 
of the nonlinear supersymmetric cr-model |j, Il8| describ- 
ing the physics of interacting diffusive modes. The latter 
problem can be solved using the transfer-matrix tech- 
nique introduced by Efetov and Larkin ^19. |. The idea 
of the method is to reduce the evaluation of the func- 
tional integral over the superfield Q{x) to the solution of 
a differential equation, in analogy with constructing the 
Schrodinger equation from the Feynman path integral. 
However the resulting transfer-matrix equation at finite 
frequencies uj appears to be too complicated (even in the 
simplest unitary case), so that no results about correla- 
tions of different eigenfunctions had ever been obtained 
in the nonperturbative localized regime, lo < A^, where 
Aj is the level spacing at the scale of the localization 
length. Instead, one can solve the transfer-matrix equa- 
tion in the limit of w — > -\-iO thus extracting information 
about statistics of the same wave function. In this way 
one can find the localization length for various symme- 
try classes 19] and calculate the whole distribution func- 
tion of a certain eigenfunction's intensity |-0„(r)p |20j. 
Similarly, the low-frequency behavior of the dissipative 
conductivity in the QID case is not known up to now. 

Thus we see that our knowledge about QID localiza- 
tion is far from being complete. This is a challenge since 
it is quasi- rather than strictly one-dimensional localiza- 
tion, which is important for various applications. The 
main technical problem lies in the complicated form of 
the transfer-matrix equations which did not allow to get 
any results concerning information about different eigen- 
functions in the localized regime, lo < A^. 

In this Letter we make the first step towards the theory 
of QID localization at finite frequencies by calculating 
the correlation function of the LDOS 



R{uj;ri,r2) = i^ ^{Pe{ri)pe+u{r2)) 



(1) 



in a thick disordered wire in a magnetic field (unitary 
symmetry); pe(r) = J^n IV'n(i")p(5(e ~ e„), and i^ is the 
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i?™'"«(w) = 1 + A{lu) + B{uj) 




^closc(-^) = 1+A(cj) 



FIG. 1: The LDOS correlation function vs. u/ A^: at coinci- 
dent points, i?'^°'"=(tj), and at fc-^ < |ri -r2| < \, i?'='°"°((.j). 



three-dimensional (3D) density of states for spinless par- 
ticles. In the present paper we restrict ourselves to the 
limit of sufficiently small spatial separations |ri — r2| ^ ^, 
when only the zero mode of the transfer-matrix equation 
is relevant. This zero mode is found exactly using the 
mapping onto an auxiliary 3D Coulomb problem. The 
explicit expression (|14|) for the zero-mode wave function 
is the main technical achievement of this Letter. 

Our result for the LDOS correlation function in the 
limit |ri — r2| ^ ^ has the form 



which is smaller than 1 indicating level repulsion (bottom 
curve in Fig.^. However, this repulsion is weak enough, 
^ciosc^-Q-j _ 2/3^ compared to the perfect level repulsion, 
i?(0) = 0, in the random matrix theory. Surprisingly, 
this value of 2/3 exactly coincides with the result [l5j 
obtained for the strictly ID geometry for cj ^ in the 
equivalent limit kp ^ r <C ■fiD- 

We start the technical section of this Letter by repre- 
senting the correlation function (^ in terms of the re- 
tarded and advanced Green functions: 

P, , 1 , Re(G,V(ri,ri)G'^(r2,r2)) 
i?(^;ri,r2)^- + ^^^^ . (3) 

Then we write G^ and G as the integrals over super- 
fields and perform the standard sequence of steps leading 
to the QID unitary cr-model 18] formulated in terms of 
the supermatrix field Q acting in the direct product of the 
Fermi-Bose (FB) and Retarded-Advanced (RA) spaces. 
The resulting expression has the form 



i?(c.;ri,r2) = \ - ^Re J P[Q]e~'^Q^DQ{x), 
with the usual diffusive action and the preexponent: 

S[Q] = ^ str / [D{S/Q{x)f + 2iu;AQ{x)] dx, 



(4) 



(5) 



lAA/ 



-iRA/ 



-lAR/ 



i?(w;ri,r2) = 1 + A{lj) + k{ri,r2)B{Lj). (2) ^IQ] = QBBixi)QBBixi) + k{ri,r2)QBBixi)QBBi^i) 



The factor fc(ri,r2) — (ImG'"'^(ri .r-?,))'^ /(tti')'^ accounts 
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and depends on the position of the points ri and r2 with 
respect to the boundaries of the wire. fc(ri , r2) is equal to 
1 at coincident points and decays fast as |ri — r2| exceeds 
the Fermi wavelength 21^ /kp. 

The universal functions A{uj) and B{uj) determined by 
slow diffusive modes are insensitive to a particular geome- 
try of the wire. They depend on the single dimensionless 
parameter w/Aj, where A^ = (Air^Di/i)'^ is the level 
spacing within the localization region {D is the diffusion 
coefficient, i^i = Av, and A ^ Nk^"^ is the wire cross- 
section). The monotonous functions Aiuj) and -B(w) are 
given explicitly by Eq. (|15|l and (|17|l . The function A{Ld) 
is negative with A(oo) = and A{Q) = —1/3, whereas 
B{uj) is positive with B{oo) = and 5(0) = 2/3. 

At coincident points, the correlator i?™'"''((jj) = 1 + 
A{ijj)+ B{u}) is always larger than 1 signaling short-scale 
level attraction (top curve in Fig. ^). It starts with 
^coins^Q^ _ ^^3 jj^ ^j^g deeply localized regime, has a 
maximum in the crossover region u> ~ A^, and slowly (as 
y^^'K^JuJ) reaches the uncorrelated limit R{oo) = 1 in the 
metallic regime (w ^ A^). This behavior should be con- 
trasted with the w-independence of the same correlator 



Dcoins 



(w) = 1 for strictly ID chains jig . 



For sufficiently close but different points, kp ^ |ri — 
r2| «C ^, the correlation function |(2Jl is nearly indepen- 
dent of |ri - rsl and equal to i?'='°"°(cj) = 1 + A{uj) 



(6) 

where str(. . . ) is the supertrace, A = diag(l, —1) is a ma- 
trix in the RA space, and Q"^ denotes a single element of 
the Q matrix. In deriving P[Q] we have assumed that the 
field Q does not fluctuate between the points xi and X2, 
which is true provided that |ri — r2| ^ min(^,^-y/A^/ci;). 
As we will be mainly concerned with the localized regime, 
u! < A^, we will refer to this relation as |ri — r2| <C C. 

Since the preexponent © depends only on Q in the 
single point Xi, the functional integral l@J can be written 
as the integral over the single supermatrix Q = Q{xi): 
/P[Q]e-^[QlDQ(x) = /P(Q)*2(g)dQ, where *(Q) = 
Io(x )^nS~^^^^DQ{x > xi). In the transfer-matrix ap- 
proach, the function ^(Q) is obtained as the zero-energy 
solution of an appropriate supermatrix Hamiltonian [l9j . 
In Efetov's parameterization ^4, 18], *(Q) = ^(A, Ai) 
depends only on two variables: A and Ai, parameterizing 
the FF and BB sectors of the supermatrix Q, respec- 
tively. Taking the integrals over the other variables of 
the parameterization we obtain the following expressions 
for the functions defined in Eq. (01: 

^(w) = Irg /"*2^A,Ai)dAdAi, (7) 

B{Lu) = ^ReJ^-^^^{X,X,)dXdXu (8) 

where the integrals are taken over the strip — 1 < A < 1 
and Ai > 1. 
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FIG. 2: To the construction of elliptic coordinates. 



The function 'I'(A, Ai) is the zero mode {H'i' ~ 0) of 
the transfer- matrix Hamiltonian |l8l |ljj : 
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(Ai-A)^ 
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Xj-l d 



dXi (Ai - A)2 dXi 



4A, 



■(Ai-A) (9) 



with the boundary condition ^(1, 1) = 1. 

The Hamiltonian (j^ without the potential term ex 
(Ai — A) had been studied in detail |22l |23j. The zero 
mode in the deeply localized regime, cj <C Aj, had been 
obtained in Ref. [l2|. In this limit it depends only on 
the noncompact bosonic variable Ai and has the form 
^(A,Ai) « ^-2i{uj/A^)Xi Ki{^~2i{Lu/A^)Xi), with 
Ki being the MacDonald function. 

Now we show that the zero mode of the Hamilto- 
nian can be found exactly for an arbitrary uj. The 
region of variation of the parameters A and Ai suggests us 
to consider them as elliptic coordinates on the half-plane 
y > of an auxiliary plane (x, y): A = (r— ri)/2 and Ai = 
(r-fri)/2, where r = y/x'^ + y'^ and ri = a/(x — 2)^ -|- y'^^ 
see Fig. 13 Changing variables from (A, Ai) to {x,y) and 
introducing a new function <i>: ^ = ri$, we obtain 



iJ* 



dx'^ 



aj/2 



1 d 



+ 



y dy 2A4 



$ = 0. (10) 



The differential part of this operator resembles the 
Laplace operator written in cylindrical coordinates. Ro- 
tating the upper half-plane around the x axis by the angle 
(fi we span the whole 3D space (the function $ is evidently 
independent of ip). The initial condition ^(1,1) = 1 
transforms to <i>(ri -^ 0) = 1/ri, which can be easily 
incorporated into the equation for $: 



[V^ + iuj/2A^r] $(r) = -4TrS^^\r - r'), 



(11) 



where r' is the 3D vector pointing from the origin to 
the point (2,0,0), so that |r — r'| = ri. Comparing with 
the equation {V^ - 2a/r + P}Gkir,r') = 6'^^^r - r'), 
which defines the Green function in the Coulomb po- 
tential a/r at energy k^/2, we immediately conclude 
that $ = — 47rGo(r,r') is nothing but the zero-energy 
Green function in the field of an imaginary charge a — 
— i(jj/4A^. To avoid possible confusion we emphasize that 
this 'zero-energy' refers to the auxiliary Coulomb prob- 
lem, while the frequency uj can be arbitrary. We also note 
that it is only the zero mode of the Hamiltonian (O that 



can be reduced to the Coulomb problem: for a finite E, 
the equation H'i — E'^ would correspond to the motion 
in a non-central field —iuj/AA^r — E/rri. 

Surprisingly, the Coulomb Green function has a very 
simple form in the real-space representation. According 
to the seminal result by Hostler and Pratt |24| . 



Gfc(r,r') = -^ ^-Adu-d.) 

A'n:k\Y — r'| 

X W_^a/kM2{-ikti)M_,„/ks/2{~ikv), (12) 

where u = r + r' + |r — r'|, v = r4-r'— |r — r'|, and W and 
M are Whittaker's confiuent hypergeometric functions. 
Taking the fc ^ limit of Eq. ((T^ with the help of the 
asymptotic formulae for Whittaker's functions ,2^], we 
obtain the Green function at zero energy: 

r, f /^ {du - d,)^ Ki{2y/^) ^ h{2y/^) 
Go(r,r) = ,^-^^-^^ . (13) 

In order to get the function ^(A, Ai) one has to multi- 
ply $ = — 47rGo(r, r') by ri that cancels the factor |r — r'| 
in the denominator of Eq. (|13() . and express u and v in 
terms of A and Ai: u — 2(Ai + 1), t; = 2(A + 1). Thus we 
obtain the explicit form of the zero mode of the Hamil- 
tonian valid for an arbitrary frequency uj: 

^{X,Xi)^Koip)qh{q)+pK,{p)Io{q), (14) 



where we have denoted p = y/—2i{uj/A^){Xi + 1) and 
q — y/—2i{uj/A^){X + 1). The properties of the Bessel 
functions ensure proper normalization: ^(1, 1) = 1. 

The form of Eq. H14|) indicates that the initial co- 
ordinates A and Ai were somehow 'more natural' that 
the auxiliary 3D Coulomb coordinates. Nevertheless we 
think than the mapping to the Coulomb problem is im- 
portant as it reveals the high symmetry of the transfer- 
matrix Hamiltonian (1^ re lated to the 0(4) symmetry of 
the Coulomb problein |23. EtJ . 

The knowledge of the zero mode ^(A, Ai) allows us to 
calculate A{u!) and B{u!) in Eqs. Q and © analytically. 
The integration for A(w) is done in terms of indefinite 
integrals of the product of two Bessel functions. The 
integral for B{uj) is much simplified in the representation 
of the 3D Coulomb coordinates: B{u!) — 47rRe/3(w) with 



(3{u;)^ I d^rGl{r,r') = - lim lim^^^?^' (15) 



which is calculated using the exact expression (|12|l . The 
results for A{uj) and B{uj) read (k — yj —Aiuj / A^): 

A{uo) = (4/3)Re|K2[j2(^) _ /o(«:)/2(«;)] 

X [Kf{K) - Ko{K)K2iK)] - lf{K)Kf{K)Y (16) 



B{cj) ^ {4/3)Rc[h{K)Ki{K)+2l2{n)Ko{i^)], (17) 

where we have used the properties of the Bessel functions 
and the fact that k^ is imaginary {u is real) in order to 
simphfy the resuh. Exact Eqs. (|16|l and H17|l interpolate 
between the metallic regime (w » Aj), which can be 
treated perturbatively by expanding in diffusive modes, 
and the nonperturbative localized regime {cu <C A^). The 
asymptotic expressions have the form (z = lu/A^): 



A{iv) 



Biu) 



-l/3 + z^ln'^il/z}/3+..., z<l, 

-3/(64z2) + ..., z>l; 

'2/3-2z2ln(l/z)/9+..., z < 1, 

v/l/2z + ..., z>l. 



(18) 



(19) 



Now let us compare eigenfunctions' correlations in the 
ID and QID geometries. It is known that the single 
eigenf unction statistics in the two cases are closely related 
|2(]j . Namely, the statistics of the wave function envelopes 
are precisely the same while the short-scale oscillations 
are different. Having derived the LDOS correlator in the 
QID case, we are now able to compare correlations of dif- 
ferent eigenfunctions in the QID and ID cases. The most 
striking mismatch between the two problems is seen in 
the LDOS correlation function at coincident points: our 
result for R'^°^'^^{cu) is always > 1 with a nontrivial uj- 
dependence (top curve in Fig. [B, whereas Rl'^'^^{uj) = 1 



.11, V 



for the strictly ID geometry [13] • However since the 
scales shorter than the Fermi wavelength can hardly be 
resolved (and with the results being model-dependent), 
it is more instructive to compare R'^^°'^°{uj) at spatial sep- 
arations |ri — r2| larger than kp^ and smaller than the 
localization length (^ for QID, and I for ID). For the ID 
geometry, Rl^^^iuj) = 1 - (1/3) J^ dtt^e-y[{u)Tf + e] 
|l3L llq . which is definitely different from our result, 
^ciose(^-) ^ 1 + A{uj). Though i?'=i°'''=(0) = Rfg'^'iO) = 
2/3, already the next term in the small-w expansion 
[w^ ln^(A^/aj) vs. Ci;^ln(l/a;r)] demonstrates the differ- 
ence between the two problems even in the deeply local- 
ized regime. Thus we conclude that the analogy between 
ID and QID localization seen in the correlations of sin- 
gle wave functions, does not extend to the correlations of 
different wave functions. 

Finally, we mention that the same approach can be 
employed in calculating the density-density correlation 
function which determines the kinetic response of the 
wire. In the limit |ri — r2| ^ ^ we get a counterpart 
of Eq. (0) with 1 + A{uj) and -B(w) interchanged: 

(ImGf(ri,r2)ImGf+Jr2,ri)) 

= {TTvf[k{Yl,V2){l + A{lo)) + B{lo)]. (20) 

Note however that Eq. H2()|l is not sufficient to get the 
conductivity of the wire, which requires a much more 
complicated analysis of the limit |ri — r2| ~ ^ |l8lll9l|. 



To conclude, we obtained the explicit expression for the 
zero mode of the supersymmetric transfer-matrix Hamil- 
tonian of a thick disordered wire in a magnetic field. 
This allowed us to get the nonperturbative result for the 
pair correlation function of LDOS, which describes cor- 
relations between different eigenfunctions, in the limit 
|ri — Y2\ ^ C- We demonstrate that these correlations 
are different in the ID and QID cases. Generalization to 
arbitrary |ri — r2 1 and other symmetry classes will be the 
subject of further studies. 

We thank D. M. Basko, M. V. Feigel'man, I. V. Gornyi, 
D. A. Ivanov, I. V. Kolokolov, V. E. Kravtsov and A. D. 
Mirlin for fruitful discussions. The research of M. A. S. 
was supported by the Program "Quantum Macrophysics" 
of the Russian Academy of Sciences, RFBR under grant 
No. 04-02-16998, and the Russian Science Support Foun- 
dation. 



[3] 
[4] 
[5] 

[6; 

[7 

[s; 

[lo; 
[11 

[12 

[13 

[14 
[15 

[16 

[17 

[18 

[19 

[20 
[21 

[22' 
[22, 



N. F. Mott and W. D. Twose, Adv. Phys. 10, 107 (1961). 

M. E. Gertsenshtein and V. B. Vasilev, Teor. Veroy- 

atn. Primen. 4, 424 (1959); ihid. 5, 3(E) (1959) [Theor. 

Probab. Appl. 4, 391 (1959); ihid. 5, 340(E) (1959)]. 

D. J. Thouless, Phys. Rev. Lett. 39, 1167 (1977). 

K. B. Efetov, Adv. Phys. 32, 53 (1983). 

S. lida, H. A. Weidenmuller, J. A. Zuk, Ann. Phys. (NY) 

200, 219 (1990). 

Y. V. Fyodorov and A. D. Mirlin, Phys. Rev. Lett. 67, 

2405 (1991). 

G. Casati, B. V. Chirikov, J. Ford, and F. M. Izrailev, 

Lect. Notes Phys. 93, 334 (1979). 

F. M. Izrailev, Phys. Rep. 196, 299 (1990). 

S. Fishman et al, Phys. Rev. Lett. 49, 509 (1982); D. R. 

Grempel, R. E. Prange, and S. Fishman, Phys. Rev. A 

29, 1639 (1984). 

A. Altland and M. R. Zirnbauer, Phys. Rev. Lett. 77, 

4536 (1996). 

M. A. Skvortsov, D. M. Basko, V. E. Kravtsov, Pis'ma v 

ZhETF 80, 60 (2004) [JETP Lett. 80, 54 (2004)]. 

V. L. Berezinsky, Zh. Eksp. Teor. Fiz. 65, 1251 (1973) 

[Sov. Phys. JETP 38, 620 (1974)]. 

V. L. Berezinsky and L. P. Gor'kov, Zh. Eksp. Teor. Fiz. 

77, 2498 (1979) [Sov. Phys. JETP 50, 1209 (1979)]. 

N. F. Mott, Phylos. Mag. 17, 1259 (1968). 

L. P. Gor'kov, O. N. Dorokhov, and F. V. Prigara, Zh. 

Eksp. Teor. Fiz. 84, 1440 (1983) [Sov. Phys. JETP 57, 

838 (1983)]. 

U. Sivan and Y. Imry, Phys. Rev. B 35, 6074 (1987). 

O. N. Dorokhov, Pis'ma v Zh. Eksp. Teor. Fiz. 36, 259 

(1982) [Sov. Phys. JETP Lett. 36, 318 (1982)]. 

K. B. Efetov, Super-symmetry in Disorder and Chaos 

(Cambridge University Press, New York, 1997). 

K. B. Efetov and A. L Larkin, Zh. Eksp. Teor. Fiz. 85, 

764 (1983) [Sov. Phys. JETP 58, 444 (1983)]. 

A. D. Mirlin, Phys. Rep. 326, 259 (2000). 

Ya. M. Blanter and A. D. Mirlin, Phys. Rev. E 55, 6514 
(1997). 

B. Rejaei, Phys. Rev. B 53, R13235 (1996). 

A. D. Mirlin, A. Miiller-Groeling, M. R. Zirnbauer, Ann. 



Phys. (NY) 236, 325 (1994). [26] V. Fock, Z. Physik 98, 145 (1935). 

[24] L. Hostler and R. H. Pratt, Phys. Rev. Lett. 10, 469 [27] M. Bander and C. Itzykson, Rev. Mod. Phys. 38, 320 

(1963). (1966); ibid. 38, 346 (1966). 

[25] H. Bateman, A. Erdelyi, Higher transcendental functions, 

Vol. 1 (McGraw-Hill, New York, 1953). 



